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Abstract 

In this paper we derive the equations characterizing the boundary layer which describes the 
transition of the distribution function of a gas of weakly interacting bosons to the distribution 
function of the gas in the presence of a Bose-Einstein condensate. To this end, we first rederive 
the classical Uehling-Uhlenbeck equation taking as a starting point the dynamics of a system of 
many weakly interacting quantum particles. The solutions of the Uehling-Uhlenbeck equation 
yield blow-up in finite time. Near the blow-up time the approximations used to derive the 
Uehling-Uhlenbeck equation break down. We derive the set of equations that describe the building 
of correlations and the onset of quantum interference effects for the many-particle hamiltonian 
system under the assumption that the blow-up for the Uehling-Uhlenbeck equation takes place in 
a self-similar form. 
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I. INTRODUCTION. 



The kinetic equation that describes the evolution of a rarefied system of bosons was 
obtained in by L. W. Nordheim and in 26] by E. A. Uehhng and G. E. Uhlenbeck. In 



the hmit of Born's approximation, the resulting equation, also known as quantum Boltzmann 
equation is the following: 



If + -V./(pi,x,t) = C(/,/) 
at m 



(1) 



C{f,f) 



4:TTg^ f dp2 



h J (27r7i)3 J i2Tihf J {2Tihf 

X ?[/](Pl,P2;P3,P4,i) (2) 

g[/](Pl,P2;P3,P4,t) = [/(P3)/(P4)(1 + /(P1))(1 + /(P2))- 

-/(Pl)/(P2)(1+/(P3))(1 + /(P4))] (3) 

with g = ATcah'^/m where m is the mass of the particles, e{p) = |pp/2m is the energy of the 
particles and a is the s-wave scattering length (cf. 28|] for example). 

The starting point for the derivation of a kinetic system of equations is a set of equa- 
tions describing the dynamics of a system of particles. In the case of weakly interacting 
bosons such dynamics is described by the Schrodinger equation for a system of interact- 
ing particles. Assuming that the interaction between different particles are weak enough it 
is possible to obtain suitable evolution equations for the one-particle distribution function 
using a perturbative method. For classical particles, this has been rnathematically proved 
for short times (cf. [19]) or globally in time for special situations (cf. [14]). This is the stan- 



dard method used to derive the Uehling-Uhlenbeck equation (see for example the c 
monographies [l, 3], as well as the more mathematically oriented approaches in 



assical 



and 



24|). Similar arguments for the Fermionic case may be found in 13| and |12| . 



It turns out that the solutions of ([I])-(l3l) can develop singularities in finite time, as it 
has been obtained in the numerical simulations for spatially homogeneous distributions of 
particles in 22|], 23(] as well as in 18|]. The interpretation of this blow up phenomenon, given 
by the authors of these papers, is that such event corresponds to the formation of the B-E 
condensate. As we will see in this paper, the derivation of the U-U equation, taking as a 
starting point a quantum many particles system, is not valid near the time of the formation 
of the singularity. 
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On the other hand, the quantum dynamics of the particles in the condensate is described 



by the Gross Pitaevskii equation (cf. 



mm 



28j). A rigorous mathematical proof of the 
for the stationary case and in 9| 
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precise formulation of this fact has been obtained in 
for the non equilibrium case with short range interactions in suitable scaling limits. 

We are interested in the process formation of the condensate, a q uestion which has 
already been considered by several authors, see for example 6|, 16|, 25||. Our main goal is 
to describe in a detailed manner how the transition between the kinetic regime described by 
the U-U equation, and the quantum regime described by means of a nonlinear Schrodinger 
equation. The relevant non-dimensional parameter is the quotient AEAt/h, where AE is 
a characteristic value of the energy and AT is a characteristic time scale for the variation 
of the density distributions. Interference effects cannot be ignored if this non-dimensional 
number becomes of order one. We derive then the equation of the boundary layer which 
should describe in detail the transition from the distribution function at the critical time to 
the density function in which the condensate is present. 



Part of the arguments presented in this paper are standard in statistical physics. 
For example we make extensive use of the second quantification formalism in Section 11. We 
also use the classical BBGKY hierarchy for quantum particles in Section III. The arguments 
used in these two Sections are also used in the theory of weak turbulence and in general in 
the derivation of kinetic equations for weakly interacting waves (cf. ^, 27]). Part of the 

281]. The main result of 



arguments of Section IV are reminiscent of similar ones in Mi and 
the paper is contained in Sections V and VI. We have however included the results in the 
previous Sections in order to fix the notation and explain the limit under consideration. 



II. THE N-PARTICLES SYSTEM AND THE SECOND QUANTIFICATION FOR- 
MALISM . 

We start recalling the classical derivation of the U-U equation that takes as starting 
point the study of the dynamics of a quantum A^— particle system by means of the second 
quantification formalism. This will allow us to precise the assumptions in which such a 
derivation is based. 

Let us assume that we have N quantum particles contained in a box Q = ]0, LY . We will 
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denote the density of particles as = -p-. We will suppose also that the particles interact 
by means of pair potentials. The hamiltonian of the system is then given by: 

Hn = Hq^n + HiN (4) 

where: 

TV , Af TV 



The evolution of the system is given by means of Heisenberg's equation for the density 
matrix: 

ihdtp=[HN,p] (6) 

p(0)=Po. (7) 

The precise form of the initial density matrix that characterizes the initial state of the system 
will be given later. Let us precise the order of magnitude of the several parameters arising 
in the system. There are two main characteristic length scales in the problem, namely the 
average distance between particles d = and the De Broglie length that is given by 

A = ^ = ^2mkBT ' assumption pq = y/2mkBT is just an estimate of a typical particle 
momentum. The temperature T is not a true thermodynamic temperature, because the 
system is not in equilibrium, but it is just a measure of the characteristic energy for the gas 
particles. 

On the other hand, Born's approximation means in mathematical terms that we may 
assume that the interaction potential between particles V is smooth, but after deriving a 
set of kinetic equations we will take the limit V g 6 (x), where g is defined just after the 
formula ([3]). 

The main assumptions on the physical parameters that we use in this paper are the 
following: 

N » 1, L >> 1, n = — = constant (8) 
A ~ (9) 

V^~lP«i- (10) 

Assumption ([8]) is the usual Thermodynamic limit assumption that ensures that there are 
not boundary effects on the resulting equations. Assumption (|TOl) means that the particle 



interactions are weak and allows us to derive a kinetic equation for the evolution of the 
distribution of particles. Finally assumption ([9]) means that the particle densities are large 
enough to allow for the formation of B-E condensate. This is related to the fact that the 
kinetic equation obtained under the assumption can yield blow up in finite time 

We will first obtain a set of kinetic equations describing the evolution of the solutions of 
([6]), ([7]) in the limit defined by (!8|)- (fT0|) . This problem was solved by Nordheim and Uehling 
& Uhlenbeck under the implicit additional hypothesis of the boundedness for the solution 
of the resulting kinetic equation. However, such assumption fails, because the solution of 
the limit equation blows up in finite time as it has been seen in the numerical simulations 



of 



18|, 



22 



231]. Therefore the Uehling- Uhlenbeck equation is not the correct limit for the 



system of particles under consideration in the limit fl5])- ffTU]) . if the time t is close to the time 
of formation of a singularity. The main goal of this paper is to obtain a new kinetic equation 
describing the distribution of particles during the formation of the condensate. 



A. Second quantification formalism. 

In order to study the particle system in the limit fl8|)- f|T0|) we will use the formalism 
of the second quantification. Most of the computations in this Subsection are standard, but 
we will reproduce them for the reader's convenience. We will assume that the hamiltonian 
Hn in (Hj), ([5]) acts in the Hilbert space: 

TV 

Hn ^ <S) i^' ^^^T (11) 

n=0 

By definiteness we will assume that the wave functions satisfy periodic boundary con- 
ditions in Q. Homogeneous Dirichlet boundary conditions would work similarly. For peri- 
odic boundary conditions the eigenvalues of the momentum operator for a single particle 
Pk = -ihd^^ are given by: 

ij 

We will denote the normalized eigenf unctions associated to the operator -ffo,Af as: 

|...,n£,....) , £ez3 
where is the number of particles in the state £. 
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For notational convenience we will use also the following alternative way of writing these 
eigenfunctions: 

|...,n^,....) = \n{£)), where n : ^ N = {0,1,2,...} 

This notation will be convenient to write in a short manner all the possible choices of 
occupation numbers. 

We introduce the well known annihilation and creation operators a^, , Ni whose action 
on these eigenfunctions is given by: 

|...,n^, ....) = a/^I---,"'^ - 1, ••••) , ^ e (12) 
|..., ....) = + 1 ni + 1, ....) , £ G (13) 
Ni ....) = aja^ \--,ni, ....) = rii |...,n£, ....) (14) 

These operators satisfy the commutation relations: 

[afc,a|] =4/, = [al,aj] =0. (15) 

We now define the annihilation and creation operators of a particle at the point x of by 
means of: 

^[x) = -j^aie i , V (3^) = -T 2^ l (16) 

ZTo L2 — 



Notice that using ( ITSjl : 

[■0 (x) , ^+ (y)] = ^ e'""'*i^" =5{x-y) 

\^{x),^{y)\ = [ip+{x),4j+{y)] =0 
Using all these operators we can rewrite the operator Hq^n as: 

Ho,N = ^^ = Y1 ^i^i^e 

where: 



Taking the gradient of f|T6l) we obtain: whence: 

^o,7v = / VV^+ (x) VtA (x) rfx (17) 



On the other hand: 



Hi,N = ^ I dxi I dx2V {xi - X2) ip (xi) ip^ {X2) ip (xi) ip {X2) (18) 



in Jn 

We define the distribution functions : 



4^ (xi, xf yi, ...,ym) = Tr (p^+ (yi) ^+ (1/2) -ip^ (ym) (a^i) (xa) .-i) {xj)) (19) 

The computation of the evolution equations for the functions fj^m is standard (see for ex- 
ample [l|) . Using ([6]) it follows that: 

ihdtf.^m = Tr (p [^+ (yi) ^+ (2/2) (2/™) ^ (xi) ^ {X2) ..4 {xj) , i/^] ) (20) 



On the other hand, we can compute the commutator in ( 1201) to obtain the following 
evolution equation for the distribution functions fj^m '■ 



ihdtfj 



,s=l 



+ / di 



n 



J2^i^-^k)-J2^{^-yk) 



.k=l 



k=l 



fj+l,m+l {^1, Xj, yi, ■■■,ym,0 + 



1 

+ 2 



s=l , Z^T^s 



k=l s=l 



(21) 



B. On the choice of the initial data. 

In order to solve the system of equations (pT!) we must prescribe suitable initial data. We 
will assume that the initial matrix density p (0) satisfies: 



where: 



p(0)=Po = ^ J2 Po{z,e;n)\n) 



.n 



n 



00 



n] 



(22) 

(23) 
(24) 



has been chosen in order to have Tr (po) = 1, and where: 

^ ® (1^1^) • 

Choosing the initial data as in (l22!l . (123|1 we are assuming that the particles are independently 
and homogeneously distributed in space according to the distribution O (■) in the space of 
energy. Since we use a macrocanonical distribution the number of variables ia stochastic 
variable. In the thermodynamic limit the fluctuations in the number of particles can be 
expected to disappear as it is usual in statistical physics. The value of z is chosen to obtain 
a given average number of particles iV for the distribution. Therefore: 



(iV)=Tr(poiV) = i Po{z,Q-n)N{n) = z 



9(log(Q)) 



where N {n) = J^eez^ n {£). 

Instead of analyzing the original system 1^, ([7]), we will study the equivalent system of 
equations (12T!) that is more convenient to use perturbative arguments. Due to (fT9!) we must 
solve these equations with initial data: 

fj,m;0 {xi, Xjiyi, ...,ym) = (25) 
Tr (po^^ (xi) ip-^ (xs) ...ip^ (xj) ijj (yi) ^/j (1/2) •••^ (l/m)) • 

Using ([22D, (123]) we obtain: 

fj,m;o{xi,...,Xj;yi,...,ym)=0 if jj^m (26) 

The evolution equations (I2T1) . with initial data (1261) admit a solution satisfying: 

fj^mixu...,Xj;yu...,ym;t) = if j^m 

Therefore, we can restrict our study to the functions: 

Fk (xi, ...,Xk; yi, ...,yk; t) = fk,k {xi, Xk, yi, ...,yk; t) 

On the other hand we can compute the initial distribution 

Fk,o (xi, Xk] yi, yk) = Fk (xi, x^; yi, yk] 0) that due to and ((251) is given by: 
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Fk,o{xi, ...,Xk;yi, ...,yk) 



I3k 

Tr {poai^ai^...ai^aj,...aj^) . 



(27) 



By ( 122]) we then have: 



We now compute the terms (n| a^^al^. . .a^^aj^. . .aj^^ \n) . We will assume for the moment that 
all the terms ji, Ja; are different. In this case, the matrix element (n| a^_^aj'^. . .a^^aj-^. . .aj^ \n) 
is not zero if and only if the coefficients ii, ...,ik are the same as the ji, or a permutation 
of them. Therefore: 

k 



{n\ a^^al^. . .a^^aj^. . .ajj_ \n) 



s=l 



]a(2)/2---^ja(k),t 



(28) 



where S"^ is the group of permutations of the elements {1, ...,k} . If some of the coefficients 
ji, ...jjk are repeated, the scalar product on the left hand side of (128!) can be bounded, using 



as 



. Therefore, using 



standard statistical physics computations (cf. c 
formula (!28l) again and replacing some sums by Riemann integrals, we can approximate (!27|) 
in the limit L ^ oo as 

k 



Ffc^o {xi, ...,Xk; yi, ...,yk) = ^ 11^° (^"^ ~ ^'^M' ^) 

o-g5fe m=l 



(29) 



where: 



Foiy;z) 



zQ (e(O)e^"^^ 

1 - ^0 (^(0) , 



(30) 
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III. THE SMALL CORRELATIONS APPROXIMATION. 



A. Non dimensional equations. 



Summarizing, we have reduced the problem to the following system of equations (cf. fl?Il) . 



ihdtFk (xi, Xk, yi, Vk, t) = Ai + A2 + A3; k = 1,2, 



^1 = ( Yl [^^^ ~ ^y^^ ) (^1' ^) 

Ffc+i (xi, Xfc, ^; Vk, ^; t) (33) 



2 = I 

n 



k 



(31) 
(32) 



Fk{xi,...,Xk;yi,...,yk;t) (34) 



(35) 



^[V(e-x,)-l^(e-%-)] 

^3 = ^ E E [V^(^.-:r.)-V^ (%•-!/.)] 
Li=i s=i , jf^s 

with initial data: 

k 

Fk,o {xi, Xk] yi, ...,yk) = E 11 "^o ~ a;^(m); ^) • 

o-eS* rn=l 

The two first equations of this hierarchy are: 

ihdtFi {xi;yi;t) = -— (A^^ - Ay J Fi {xi, yi, t) 
2m 

+ [ d^[Vi^- X,) -V{^- y,)] [F, (xi; y,; t) F, (C; t) + F, {x,; t) F, (^; y,; t)] 
Jn 

+ [ d^[V{^~x,)-Vi^-y^)]F2{x^,C,y^,C,t) (36) 
Jn 



ihdtF2 {xi,X2; yi,y2; t) = -— (A^^ + A^^ - ^yi - ^2) F2 (xi, X2; yi, 2/2; t) + 
+ di[V{i- x,yv (e - %•)] F3 (xi, X2, 1/2, e; t) + 

+ [V (xi - X2) - F (yi - y2)] F2 (xi, X2; 1/1,1/2; t) (37) 

In order to understand more clearly the limit that we are considering, we introduce the 
non-dimensional variables: 

x = \x, Vix) = ^Vix), Fkix)=i^j Fk{x), t=-^i, p=-p (38) 



9 = e 



:d' 



(39) 
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where, due to ffTOj) . e is a small parameter and by assumption, the potential V{x) is now of 
order one. Our choice of time scale is due to the fact that we want to obtain, in the limit 
e ^ 0, an equation in which the particle density varies in times t of order one. Notice that 
{x,p) — > {x,p) is not a canonical transformation but transforms a "quantum cell" in the 
phase space of volume h into another cell of volume one. 

For the sake of simplicity we drop the hats in all the variables x, p, i, and V as 
it is customary in the computations of asymptotic expansions. The system ( 136|) - (!37|) then 
becomes 

idtFi (xi; yi; t) = (A^^ - A^jFi (xi; yi, t) 

+ - [ [V (e -Xi)-V{i- y,)] [Fi (xi; y,- 1) {i; ^; t) + F, (xi; t) F, i^; y,; t)] 

+ - [ d^[V{^-x,)~V{^-y,)]F,{x,,^;y,,^;t) 

(40) 

idtF2 (xi, X2; yi, 1/2; t) = (A^^ + A^^ - A^^ - Aj^J F2 (xi, xs; yi, 2/2; t) + 
2 

+-y, [V (e - x,)-V (e - y,)] F3 (xi, X2, 1/1, 1/2, t) + 
' ,=1 



1 fd^ ^ 



[y{xi-X2)-V{yi-y2)]F2{xuX2;yuy2;t) (41) 



B. Small correlations limit. 

Our aim is to obtain closure relations for the functions F^ by means of a perturbative 
argument. 

Notice that in the absence of potential the system of equations ( 1401) . (HTl) might be 
explicitly solved and the resulting solutions have the form: 

k 

Ffc (xi, ...,Xfc;|/i, ...,|/fc;t) = ^ JJ Fi (x^(„), t) (42) 

Moreover, the function Fq in fl30p is invariant under spatial translations whence 
Fo(xi;?/i) = -Fq (xi —yi). Since the system of equations (!3T]) - (!MI) are also invariant un- 
der spatial translations it follows that Fi (xi; yi, t) = Fi (xi — yi] t) for any t > 0. 
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Notice that in this case we can think in the solutions of this form as "uncorrelated" 
solutions, although in a strict mathematical sense the corresponding probability distributions 
are not uncorrelated, but the only correlations between particles whose distribution is given 
by fH21) would be the ones due to the symmetry of the wave functions due to the bosonic 
character of the particles (cf. the discussion in pj]). In any case the approximation (1421) is a 
convenient starting point for the computation of the solutions of (l3T]) - (l35|l in a perturbative 
manner. We define the correlation functions Gk by means of the identity: 

Fk (xi, rcfc; 1/1, Vk] t) = Gk (xi, Xk, yi, Vk; t) + 
+ Fk{xi,...,Xk;yi,...,yk;t). (43) 

where we have defined 

k 

Fk (xi, ...,Xk;yu ■■■,yk;t)=^ JJFi (x^(„), y^; t) . (44) 

It is possible to derive a kinetic approximation for (136!) (137|) under the following small 
correlation assumptions 



\Gk\ « n l^il- (45) 



m=l 



Indeed, under this assumption we obtain, plugging P3l) into (I36l) - fl37j) : 

idtFi {xi; yi; t) = (A^.^ - A^jFi (xi; yi; t) 

+ - [ [V (e - xi) - F (e - yi)] [Fi (xi; yi; t) F, ^- 1) + Fi (xi; t) Fi (e; yv, t)] 

+ - / di[V{i-x,)-V{i-yr)]F^{x,,i-y,,i-t) 
f Jo. 



(46) 



idtF2 (xi, X2; 1/1, 1/2; t) = (A^.^ + A^.2 - Ay^ - A^^J F2 (xi, X2; 1/1, 1/2; t) + 



1 /■ 

+- V / cie [V^ (e - (^ - y^)] F3 (xi, X2, Z/i, 2/2, + 

' ,=1 



1 ^ 



+-(^) [^(a;i -X2) -K(yi -y2)]F2(xi,X2;2/i,2/2;t) (47) 
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The relative strength of the terms yielding correlations is of order e. This explains why in 
equation (H7|) we have approximated F2 and F3 by F2 and F3 respectively. Notice that we 
have kept all the terms in the equation (1461) and only terms of order 1/e or larger in fl47l) . 

We now compute the evolution equation for ^2 (xi, X2; yi, 1/2; ^) using (1^51) and the ap- 
proximation fH^ - fH7|) : 

idtG2 (xi, X2; yi, y2] t) = --^ (^^-i + ^^-2 - \i - ^2) G2 (a^i, xa; 1/1, 1/2; t) + 

+ - [t] [V (xi -X2)-V (2/1 - 2/2)] F2 (xi, X2; 2/1, 2/2; + 1) (48) 

2 

m-, = E - ^:') - ^ - %■)] ^3 (xi, X2, 1/1, 

- \y{^-xi)-y{^-y<T{i))]F2{xi,i]y^(^i),i]t)Fi {x2]ya{2)]t) 
(7e52 

- E Fi (xi; 2/^(1); t) [l^ (^ - X2) - - 2/<x(2))] i^2 (x2, l/<x(2), t) (49) 

After some computations, neglecting terms of order 0{t^) it follows that: 

H{i] t) = V{i- xi) Fi (X2; e; t) F2 (xi, 2/1, 1/2; t)-V{i- 2/1) Fi (e; 2/2; t) F2 (xi, X2; 2/1, i\ t) 
+ V{i- X2) Fi (xi; e; t) F2 (X2, 1/1, Z/2; t) - (e - 1/2) Fi (e; 2/1; t) ^2 (xi, X2; 2/2, e; i) • (50) 

Plugging (l50ll into (l48l) we obtain: 

Z(9iG2 (xi, X2; 2/1, 2/2; i) = (A^i + - Aj^j - Ay^) G2 (xi, X2; 2/1, 1/2; i) (51) 



+ - 
e 

1 

+ - 

e 



J y ^ (^ - a;i) Fi (X2; ^; t) F2 (xi, ^; 2/1, 2/2; i) - '^^ (^ - Z/i) i^i (^; I/2; t) i^2 (xi, X2; 2/1, ^; t) c?^ 



(e - X2) Fi (xi; t) F2 (X2, e; 1/1, 1/2; t)-V{^- 2/2) Fi (e; 2/1; t) F2 (xi, X2; 2/2, t) 



1 /d^ ^ 



+ [Vixi-X2)-V (2/1 - 2/2)] F2 (xi, X2; 2/1, 2/2; i) 



where, due to ([9]), (i/A is of order one. 

The equations ( H6l) and ( ISTl) provide the evolution equations for the functions Fi and G2. 
We have not used yet Born's approximation which in the variables that we are using reduces 
to 

V{x) = 6 (x) (52) 
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Using (152!) we obtain: 



idtFi {xi;yi;t) = — - (A^., - Aj, J Fi (xi; yi]t) + 



+ -g [G2 ixi,xi; yi, xi, t) - G2 (xi, yi; 2/1, yi, t)] (53) 



idtG2 {xi,X2] yi, 2/2; t) = — - (A^^ + A^^ - A^^ - A^J G2 (xi, X2; yi, 2/2; t) 



1 

e 

1 

+ -^? 
e 



Fi (x2; xi; t) F2 (xi, xi; yi, 7/2; t) - Fi {yi, y2] t) F2 (xi, X2; yi, 2/1; t) 



+ 



Fi (xi; X2; t) F2 (X2, X2; 1/1,1/2; ^) - Fi (1/2; 2/1; t) i^2 (a^i, 2^2; I/2, 2/2; i) 

+-f^j 5- - 3^2) - 5(2/1 - 1/2)] i^2 0:2; 1/2, 1/2; (54) 

The invariance of the initial distribution Fq (xi; yi) under spatial translations imply that, 
with a slight abuse of language, the solutions of ( l53l) . ( l54l) have the form: 



Fi (xi;?/i;)f:) = Fi (xi - yi]t) 
G2 (xi, X2; 2/1, 2/2; t) = G2 {xi - yi, X2 - 2/1; 0, y2 - 2/1; t) 



(55) 
(56) 



Under these assumptions the equations ( l53l) . (15^ reduce to: 



i^tFi (xi - yi; t) = [G2 (xi, xi; yi, xi;t) - G2 (xi, yi, 2/1,2/1; ^)] (57) 

1 1 

idtG2 (xi, X2; 2/1,2/2; i) = — ^ (Aa;i + - Ay^ - Aj^J G2 (xi, X2; 2/1, 2/2; t) + - Q [Fi] 



(58) 



where: 



Q [Fi] (xi, X2; 2/1, 2/2; t) = Fi (x2; xi; t) F2 (xi, Xi; yi, ?/2; t) - Fi {yi; 7/2; t) F2 (xi, X2; 2/1, yv, t) 



+ 



Due to (135!) we have: 



Fi (xi; X2; t) F2 (x2, X2; 2/1,2/2; i) - i^i (2/2; 2/1; if:) ^2 (xi, X2; 2/2, 2/2; i) 

+ (^) [.(xi-X2)-.(2/i-2/2)]F2(xi,X2;2/2,2/2;t) 

(59) 

G'2(a;i,X2;2/i,2/2;0) =0. (60) 
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The system of equations (l571) - (l60ll will be our starting point for the description of the 
Bose gas in which we are interested. Notice that it is a closed system of partial differential 
equations. 



IV. THE PROBLEM IN THE PHASE SPACE. 



The function that describes the one-particle density in the phase space in quantum prob- 
lems is the Wigner transform of Fi {xi; yi;t) . Such a function is defined as: 

A {x, p; t) = -i- / Fi (a; + C; X - C; t) e'^^dC (61) 

where the normalization constant in (|6T1) has been chosen in order to have 

j fi (x, p) dxdp = N. 

In the spatially homogeneous case we have Fi {x + x — t) = Fi {(; t) due to fl56|) . There- 
fore, (l6Tll reduces to the Fourier transform: 

A {x, p; t) = h {p, t) = -^f F, (C; t) e^PrfC (62) 

In order to obtain the evolution equation for fi {p, t) we then take the Fourier transform 

of dssD, dsni): 

tdth (P, t) = -V / (C, C; 0, C; t) - (C, O; o, O; t)] e^^dC (63) 

On the other hand we have the following Fourier representation for the function 

G2 ixi,x2;yi,y2;t): 

92 (^1, 6; Vi,V2; t) = ! dxidx2dy^dy2e'^^'^'+^^-^-^'y^-^^y^^G2 (xi, 2/1, 1/2; t) (64) 

(27r) J(R3)4 

G2 (xi, X2; yi, y2, t)= [ d^^d^2dmdv2e-'^^''''+^'"'y'-'''y''^g2 (6,6; Vi, V2; t) . (65) 
Let us write: 

w (6, 6; Vu V2; t) = / dxidx2dyidy2e'^^'^^+^'^'-^'y'-^'y'^Q [Fi] (xi, X2; yu 2/2; t) 

(27r) J(R3f 

(66) 

Taking the Fourier transform of (l58l) and using (165|) in (163!) we obtain the following system 
of equations for fi and g2 ■ 
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idtfi{p,t) = - [ d^id^2dr]idr]2[6{p-r]i) -6{p- ^i)]g2 {^1,^2; Vi,V2]t) (67) 
e J(R3)4 

-19192 (6,6; vi, V2] ^) = ^ (6) + £^ (6) - ^ ivi) - ^ im)] 92 (6, 6; ^71, V2; t) + 

+^w(6,6;^i,^2;t) (68) 

where w{^i,^2]Vi^V2]t) is as in ( l59l) . ( l66l) and the energy e(p) in the non-dimensional 
variables is e{p) = p^. To obtain a closed system for /i,5'2 it only remains to compute 
w (6, 6; ^2; ^) in terms of /i. To this end notice that fl62|) yields: 

Fi(C;t) = 1 h{p,t)e-'P<dp (69) 

Using fl55l) in the formula of Q [i^i] and plugging the final expression in fl66l) we obtain, after 
some computations: 



wi^i, 6; ^1, V2; t) = 25(6 + 6 - - ^2)g[/](6, 6; vu V2, t) 



?[/i](6,6;??i,??2,i) 



-A(6)/(6)|(^ +A(..) U^) +A(.2) 



(70) 



where we have dropped the time dependence of the function /i in the right hand side of 
rOj) . The solution 92 to (168|) is then: 

52(6, 6; r/i, t) = --5(6 + 6 - r^i - r/2) r e-^[^(«^)+^(«^)-(''^)-(''^)](*'^) x 

^ Jo 

x?[/i](6,6;m,^2;s)c/s. (71) 



where we have used that (72 (■,■;■,■; 0) = due to (160|1 and (l64ll . The Dirac measure in (17T|) 
may be simplified if we define, 

^2(6,6;^i,^2;t) = 5(6 + 6 -^1 -'?2)v5(6,6;^i,^2;t), 

from where, fl7ip gives 

¥'(6,6;r/i,r/2;t) = -- re-*[^(«^)+^(^^)-^^^^)-^(''^^^(*-^)g[/i](6,6;r]i,r/2;s)rfs (72) 

^ Jo 



16 



Plugging (1721) in (1671) and using the symmetry of q[fi] with respect to its arguments ^i, ^2, ''7i 
and ri2 we finally obtain the following equation for /i: 



dtfi {,Pi,t) 



^ -'0 



ds 



dp2dp3dp4 < cos 



ie{pi) + e{p2) - eips) - e{pi)){t - s) 



xSiPi + P2 - P3 - Pa) q[fi]iPi,P2;P3,PA; sl73) 



Non-Markovian 
for example 



Boltzmann equations have been found in several physical situations (cf. 
171], 281] and references therein). 



V. THE KINETIC LIMIT: THE UEHLING UHLENBECK EQUATION. 
A. Derivation of the Uehling Uhlenbeck equation. 

The formal derivation of the U-U equation would then proceed as follows. If we suppose 
that 



1 



'e{p) + e{p2) - e{p3) - 5(^4)) (r - a) » 1 



(74) 



a simple formal argument gives: 



cos 



ie{p) + e{p2) - e{p3) - e{pij){T - a) 



n6{T-(r)6{e{p)+e{p2) - e^ps) - e{pi)) 



(75) 



in the sense of measures, where in that formula p2 = Ps + P4 — Pi- We finally end up with 
the U-U equation: 



9tfi {pi,t) = An dp2dpsdp4S{e{pi) + e{p2) - e{pz) - ^^ipi)) 

X(5(pi + P2 - PZ - Pi) q[fl]{Pl,P2]Pz,PA]t) (76) 

Notice however that this approximation requires the condition ( 17^ . which, using the 
original physical variables, can be formulated as the quasiclassical condition AEAt » h. 
Equation (17^ is just the equation ([I])-® written in a different system of units. 
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B. The loss of validity of the kinetic approximation. 



According to the blow up scenario of Semikoz & Tkachev (cf. [22, l23|) and Pomeau et 
al. (cf. [l^, 3]) the blow up for the equation fl751) takes place in a self similar manner and 
the distribution of particles has relevant variations in the regions of the space of momentum 
p whose size rescales like the power (T — t)^ for some positive p. In order to describe this 
region by means of an equation free of parameters we look for self similar solutions of fl76l) . 
Such solutions have the form 

f{t,p) = (T - $(0, e = jjA^ (77) 



where the numerically computed value of (3 is such that (3 = 1.069 (see 18|). The function 
$, that is of order one, solves then an integro differential equation free of parameters. 
Notice also that the functional form ( l77l) tells us immediately the time scales for which the 
interference effects in f l73l) cannot be ignored, or more precisely, in dimensional variables, 
the time scale where AE At ~ h. This happens if p^(T — t) e"^ or equivalently if 

(T-t)~^. (78) 

Since p ^ {T — tY in the self similar region, we obtain that the interference effects appear 
at times 

(T - t) ~ eWT (79) 

For this time scales we have to introduce a boundary layer in order to take the interference 
effects into account. 



C. The correlations become of order one in the boundary layer. 



It turns out that in the same time scale where fj75|) starts failing, the small correlation 
approximation condition (145!) ceases being valid. Indeed, assuming the self similar behaviour 
fl77|) we obtain: 

Fi(C,t) = (27r)3 f f,{p,t)e-''^Pdp (80) 

= (27r)3(r - t)^-i/2 J $(Z)e-*^(^-*)'^rfZ (81) 
= {T -tf-^/^^{Z{T -tf). (82) 
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On the other hand, (158!) and (159|1 yield that for the boundary layer time scale: 

1 



G ~ -Ft{T - t) 



from where, we obtain, using ( 1791) : 



G ~ (T - t) 



2/3-1 



(83) 



A similar argument shows that IG^I ~ nJ^^]^|Fi| for A; > 1. It then follows that the ap- 
proximation of the system (15B]) (P7|) by system (H^ (H7|) breaks down at the time scale 

dzi. 



VI. THE BOUNDARY LAYER: ANALYTIC DESCRIPTION. 

In this Section we derive the set of equations describing the boundary layer where the 
kinetic approximation is lost. Since, as we have seen, the correlations become of order one 
in that region, we need to keep a major part of the equations in system ( 13T1) . Using the 
non-dimensional variables (138|) . that system becomes 



idtFk (xi, Xk, yi, Vk, t) = Ai + A2 + A3; /c = 1, 2, 
= \ y [ [A^-, - \s] Fk (xi, Xk] yi, ...,yk; t) 



A2 = - 

e Jn 



k 

Y^[V{^-x,)-V{^-y,)] 



Fk+i {xi, ...,Xk,^]yi, ...,yk,^]t) 



A, 



1 fd 



2e V A 



k k 

^ ^ [V{x,-x,)-V{y,-y,)] 

_j = l s = l , fe^s 



X 



The rescaling (1791) suggest to define new variables as follows: 

2/3-1 ^ 

Fk{xi, ...,Xk;yi, ...,yk;t) = e'^+i Hk{Xi, ...,Xk;Yi, ...,Yk;T) 

2 2/3 2/3 

T - t = -£2/3+1 r, = e 2/3+1 Xj, p = e2/3+i P 



(84) 
(85) 

(86) 



X Fk{xi, ...,Xk;yi, ...,yk;t) (87) 



(88) 
(89) 



Neglecting lower order terms in e and using that V{x) = S{x) we obtain that the functions 
Hk satisfy at leading order the following system: 
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^9,iJfc(Xi,...,Xfe;yi,...,n,r) = A + fc = 1,2, 



—Hk+i (Xi, Xfc, K,-; Fi, Yfc, Yf, r)] 
This system must be solved with the asymptotic condition: 

i7i(X, r, r) ~ (-r)^-^/2^((X - r)(-r)'') as r ^ -00 



(90) 



(91) 



i7fc(Xi,--- ,Xfc;Fi,--- ,n;r) ~ ^ J] i/i (X,(„),r„;r) as r ^ -00 (92) 

^g^fc m=l 

Notice that formula fl90p implies that all the correlation function Gk defined in fl43p became 
of order the order of n^^]^|Fi| in the time scale (1791) . 

This problem may be also expressed in the phase space using the Wigner transform that are 
defined as: 

^fc(Xi, ■ ■ ■ , Xfc; Pi, ■ ■ ■ , P,; r) = J dCi ■ ■ ■ dCk ^'-^^^^ x 



Plugging this into the system (19011 we obtain: 

k k 



dip. 



X |^v^fc+i(-^i5 ■ ■ ■ 5 -^fc, — 0; -Pi) • • • ) -Pj) ■ ■ ■ ) -ffc, -Pfc+i; T 

— V'fc+il-^i, ■ ■ ■ , -^fc, Xj + Cj-; Pi, . . . , Pj, ■ ■ ■ , Pfc, Pfc+i; r 
The asymptotic data as r — > — 00 are now determined by: 
^i(X; P; r) = ^i(P; r) ~ (-r)-'^-^^^ 



A;= 1,2, 



(93) 



(94) 



^,(Xi, ■ ■ ■ , Xfc; Pi, . . . , P,; r) ~ -i— ^ j dCi ■ ■ ■ dCk e' ^^^^ 



X 



X 



m=l 



(95) 
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Both systems, (!90i)- (p2l) and (p3|) - (p5|l are rather comphcated objects to study that we do 
not consider in detail in this paper. However, the solution of this problem should provide a 
clear description on how the transition from the kinetic regime to the quantum dominated 
and highly correlated regime takes place. 

It is interesting to compute the time for which the correlations appear in physical variables. 
Using (1791) . we obtain that such scale is given by: 



where T* is the time at which the Uehling Uhlenbeck equation blows up in the original 
physical units. The range of physical moments that would be described by the systems 
above ( "in the boundary layer" ) is 



2/3 



Finally the correlation lengths in this boundary layer is 



X ~ A 



a A^ 



2/3 



^3 \ 2/3+1 
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